New insights on the high-pressure phase diagram of molecular CO2 
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We report the discovery of a new molecular phase of carbon dioxide at high-pressure and high- 
temperature. Using x-ray diffraction, we identify this phase as the theoretically predicted high- 
temperature Cmca phase [Bonev et al, Phys. Rev. Lett., 91, 06550f (2003)]. Its relation with 
phase III, on one hand, and its relative stability with respect to phase IV, on the other hand, are 
discussed based on spectroscopic and melting data. The existence of this strictly molecular phase 
challenges the interpretation of phases IV and II as intermediate phases between the molecular and 
covalent-bonded forms of CO2. 



In the last years, the intensive study of simple molec- 
ular systems (H 2 , N 2 , CO, CO2 etc.) at very high pres- 
sures (P) and temperatures (T) have led to a wealth 
of remarkable discoveries, including new solid structures, 
insulator-metal transitions, symmetrization of H-bonds, 
polymerization or the formation of new covalent bonds. 
These findings have had important repercussions both 
on the fundamental level, since their understanding is 
usually challenging, and from a technological viewpoint, 
since they may open, for example, new routes for energy 
storage. 

In this context, the discovery in 1999 by Iota et al. Q 
of a non-molecular, quartz-like, phase of carbon dioxide 
by heating above 50 GPa has driven a large attention 
onto the behavior of this compound at elevated pres- 
sures. This phase was described as an extended network 
of four-fold coordinated carbon atoms, with a unusu- 
ally large bulk modulus. Subsequently, new molecular 
phases were uncovered, named II and IV (Fig. 0,0. 
From the large values of the C=0 bond lengths found for 
these two phases, and the bent geometry of the molecule 
in phase IV, they were presented as precursors to the 
non-molecular phase. This interpretation, however, is 
not supported by a theoretical study by Bonev and co- 
authors Q who found, first, that the equilibrium C=0 
bond length in any of the proposed structures is compa- 
rable to the free molecule value; and second, the stable 
structure in the domain of phase IV is not Pbcn with 
bent molecules, but a molecular Cmca structure with 
linear molecules. 

This controversy enlightens a peculiarity of carbon 
dioxide compared to other simple molecular systems, that 
is, the presence of large metastabilities when going from 
one solid phase to the other. For example, the transition 
between phase I and III is a very sluggish one, spreading 
over 10 GPa at room temperature 1 10, and sccm- 
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FIG. 1: (Color online) Experimental (black dashed lines) 
and theoretical Q (red solid lines) phase diagram of carbon 
dioxide. The theoretically predicted stable phases are indi- 
cated in italics, and the proposed experimental ones, in up- 
right characters. The dotted lines are speculated kinetic lines. 



ingly taking place through an intermediate structure • 
These metastabilities make the task of determining the 
phase lines with accuracy a difficult one and more impor- 
tantly, raise the question of which structure is actually 
the stable one at high P-T. 

In this Letter, we report a new investigation of the 
phase diagram in the pressure region below 20 GPa and 
for temperatures up to 950 K. It reveals the presence of a 
previously unobserved phase (noted CO2-VII) above 640 
K, intermediate between phases I and IV. Using x-ray 
diffraction, we determine the structure of this phase and 
show that it is the theoretically predicted high P-T Cmca 
phase. The distinction between CO2-VII and CO2-III 
is revealed by comparing their Raman spectrum. We 
also report measurements of the melting lines of this new 
phase and phase IV which allow us to determine their 
relative stability. Finally we show that the discovery of 
this phase has important consequences on the nature of 
phases IV and II. 

The experiments reported here were done in resistively 
heated diamond anvil cells. The techniques are identical 
to the ones presented in details in Ref. 
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FIG. 2: (Color online) Presently determined phase transition 
points. The melting curve of phase I reported in Ref. is 
shown as the dot-dashed line; o: melting points of phase I 
(present work); + and •: I-IV transition measured respec- 
tively on loading and unloading; A: VII to IV transition; 
□ : melting points of phase VII; ♦: melting points of phase 
IV. The experimental uncertainties on P and T are within 
the symbol sizes. The curves are fit to the data. The melt- 
ing points ofphase IV and VII are well fitted by the Simon- 
Glatzel law P : T = T [l + (P - P )/a] 1/c with T = 805 K, 
P = 11.15 GPa, a=46(8) GPa and c=3.8(6). The photo- 
graph shows the sample of phase VII (with pressure sensors) 
viewed through the diamond anvils at 770 K and 11.8 GPa. 



We first investigated the I-IV transition line between 
500 and 790 K by visual observation and Raman spec- 
troscopy, tuning the pressure around the transition in 
order to observe it both on loading and unloading. As 
Fig. [21 shows, the unloading transition points lie on a 
well defined nearly vertical line, whereas the loading 
ones are quite scattered. This hysteresis decreases with 
temperature and nearly disappears above 700 K. From 
the unloading line and the previously reported melt- 
ing curve [8j, we localize the I-IV-Fluid triple point at 
800(5) K and 11.2(1) GPa. 

For T > 640 K and P > 12 GPa, compression ofphase 
I systematically produced a new phase, noted hereafter 
CO2-VII, which in turn transformed to phase IV by fur- 
ther compression. CO2-VII was easily identifiable both 
visually and from its Raman spectrum. The transition 
was usually sharp and the new phase visibly birefrin- 
gent (Fig. |2J). Fig.|3|reports the Raman spectra collected 
along an isotherm at 719 K, following the transition from 
phase I to IV through phase VII. The characteristics of 
the Raman spectrum of CO2-VII will be analyzed below. 
At 640 K the I- VII and VII-IV transitions took place 
at 12.3 GPa and 12.55 GPa respectively. Increasing the 
temperature, the pressure domain in which phase VII is 
observed grows larger: at 726 K the two transitions were 
localized at 12.1 and 14.2 GPa respectively. Decompres- 
sion of phase VII at 705 K led to phase I at 11.8 GPa, 
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FIG. 3: (Color online) Raman spectra collected on quenched 
phase VII at 300 K(top) and along an isotherm at 719 K going 
from phase I to IV through VII (bottom). 



i.e. slightly above the IV to I transition pressure (11.4 
GPa). 

To investigate the structure of phase VII, we collected 
x-ray diffraction spectra following an isotherm at 726 K, 
using the monochromatic (A = 0.3738 A) angular dis- 
persive setup on beamline ID 2 7 of the European Syn- 
chrotron Radiation Facility in Grenoble. The images 
were collected during rotation of the cell around the uj 
axis. The observed d-spacings match very well an or- 
thorhombic unit cell with 4 molecules and the following 
lattice parameters at 12.1 GPa, 726 K: a = 4.313(1) A, 
b = 4.746(1) A, c = 5.948(1) A. All the peaks could be in- 
dexed by considering three distinct diffracting crystallites 
rotated around the c axis with respect to each other. Visi- 
ble peaks obey the following conditions: hkl : h+k = 2n; 
hkO : h,k = 2n; Okl : k.l — 2n, from which we deduce 
the space group Ccmb. This is the same group as Cmca 
but with a different choice of axes, i.e. with a and b axes 
inverted. 

The existence of a Cmca polymorph in high-pressure 
CO2 was first predicted by Kuchta and Etters using lat- 
tice dynamics calculations |l0j| . In their structure, the 
axes are in the order a < b < c and the molecules lie on 
the (be) plane, tilted by <p = 52° with respect to the c 
axis. More recent calculations based on density func- 
tional theory (DFT) [1 Q2 found that the optimized 
Cmca structure, stable at high P — T, was such that 
a > b and ip — 54°. Since we find the same axis or- 
der, we used this structure as a starting point for refine- 
ment. The latter was done against the measured inten- 
sities of 16 independent reflections originating from the 
same crystallite 0|. Convergence was obtained down 
to an i?-factor of 0.069, yielding for atomic positions, C 
on (0,0,0) and O on [0,0.212(3),0.109(1)], with isotropic 
temperature factors U of 0.0637 A 2 (C) and 0.0314 A 2 
(O). The refined structure, shown on Fig. 01 presents a 
ip angle of 54.5 (1)° , in excellent agreement with the DFT 
calculations H1I3|- The C=0 bond length is 1.13(1) A, 
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FIG. 4: (Color online) Structure of the Cmca phase VII 
at 726 K and 12.1 GPa, projected along the a axis . The 
molecules lie on the (be) plane, tilted by ip = 54.5° with re- 
spect to the c axis; in the bottom, the volume per molecule 
vs. pressure is shown along the isotherm at 726 K for phase 
I (+), VII (o) and IV (□). The solid line is the theoretically 
computed equation of state for the Cmca phase at K Q. 
The volume discontinuity is 1.2% between phase I and VII 
and 0.66% between phase VII and IV. 



i.e. slightly shorter than the one found [13j for phase I at 
1 GPa [1.1486(9) A], and again in good agreement with 
calculations. We can thus conclude that phase VII is 
the Cmca phase predicted by Bonev et al. |4J to be sta- 
ble in this P — T range. This is further corroborated by 
the excellent agreement observed between thepresently 
determined and calculated equation of state [4J (Fig. 0] 
bottom). 

It is worth pointing out that the Cmca space group is 
also the one designated for the room-temperature high- 
pressure phase III. The structure of this phase was iden- 
tified by Aoki et al. ^lj, an d later by Yoo et al. ||, as 
the one predicted by Kuchta and Etters 01- However 
we note that collecting good diffraction data on phase 
III has been so far challenging due to the highly strained 
nature of this phase. As a result, its structure is still 
not well constrained. In their work, Bonev et al. 0, 
like previous authors, assimilated the theoretical Cmca 
phase with phase III, only relocating its stability domain 
to that of phase IV (Fig. (the Cmca structure was 
found unstable at 300 K, so phase III was presented as 
metastable at this temperature) . This raises the question 
of whether phase VII is identical to phase III. First of all, 
we note that the experimental existence domains of the 
two phases are not connected in the P — T space : phase 
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FIG. 5: (Color online) Room temperature pressure evolution 
of external and internal modes in phase VII (•) and III (o). 
The mode i comes from the « 1% of 13 C 16 2 molecules. 



Ill transits to II at « 450 K whereas CO2-VII is only 
observed above 640 K. To go further, we have compared 
their Raman spectra collected during decompression at 
300 K. Phase VII could be quenched to room tempera- 
ture, which is a common property of all the high P-T 
phases of carbon dioxide. Phase III was produced by 
compression of phase I at 300 K and annealed to 440 K 
to release the stress. The pressure evolution of the Ra- 
man modes is shown in Fig. The four lattice phonon 
modes and two components of v + and v~ |l7| . observed 
in both cases, match group theory predictions for the 
Cmca structure. A « 5 cm^ 1 frequency difference be- 
tween the two solids is registered in the central lattice 
modes (b and c in the figure) and in the v + region, 
whereas the difference increases up to 10 cm -1 for the 
components of v~ . This indicates that the two structures 
are indeed different. Further work is however needed to 
clarify this question. 

By contrast with the theoretical phase diagram in 
Fig. ^ the Cmca phase is not the only one present at 
high temperatures, since it transits to phase IV at higher 
pressure. As a matter of fact, we could only obtain phase 
VII by compression of phase I and never by heating or 
decompressing phase IV. This is reminiscent of the fact 
that phase III has itself only been obtained from com- 
pression of phase I. Since this was taken as an indication 
of the metastable character of phase III, this raises the 
question of which of phase VII or IV is the stable struc- 
ture beyond the I- (IV, VII) transition line. 

To answer this, we have measured the melting curve 
of CO2 for T > 800 K, extending the determination pre- 
sented in Ref.0. A melting point was defined as the P—T 
conditions where the solid/fluid equilibrium was visually 
observed and the solid phase was identified by its Raman 
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spectrum. Melting of solid I was followed up to 808 K and 
11.7 GPa. At 810 K and 11.8 GPa the sample suddenly 
transited into pure phase VII. It was then decompressed 
in order to melt it and heated to 827 K. Increasing the 
load at this temperature produced the equilibrium be- 
tween solid VII and fluid at 12.1 GPa. We then followed 
the melting curve of phase VII up to 920 K. Compression 
of phase VII at 879 K produced phase IV at 16.8 GPa. 
Surprisingly, the latter could be decompressed down to 
the equilibrium between phase IV and the fluid, at 884 K 
and 13.2 GPa. We then measured several melting points 
of phase IV from 810 K to 946 K. As shows Fig. 2, this 
melting curve is not distinguishable from the one of phase 
VII within our experimental uncertainties (±0.2 GPa and 
±5 K at these P — T conditions) . 

The fact that the solid phase obtained from the melt 
above 827 K is phase VII confirms that it is indeed the 
stable phase at these temperatures and pressures. How- 
ever, the closeness of the melting curves shows that the 
energy difference between phase IV and VII is very small, 
since at melting the Gibbs free energy of the fluid and 
solid phases are equal. The fact that we could decom- 
press phase IV down to the solid / fluid equilibrium with- 
out reverting to phase VII gives an indication that there 
is a strong potential barrier hindering this transforma- 
tion. This also explains why we were able to measure 
the transition line between phase I and IV by compres- 
sion of phase I even in the temperature range where we 
observed phase VII: since we tuned the pressure closely 
around the transition, it is likely that part of the sample 
was still in phase IV, forcing the transition from phase I 
to IV without passing via phase VII. 

Fig. 0]shows that the VII-IV transition is accompanied 
by a small volume decrease (0.66% at 726 K). Since the 
slope of the transition line is positive, the Clapeyron's 
law imply that phase VII has a larger entropy (by about 
4 J/mole) than phase IV. This agrees with the calculation 
of Bonev et al. J4| who found that the Cmca structure is 
stabilized at high temperatures by its larger entropy due 
to the contribution of a soft acoustic phonon. Extrap- 
olating the measured melting line and VII-IV transition 
line gives a possible IV-VII-F triple point at «1030 K 
and 18.4 GPa. 

Phase VII is thus the stable solid phase beyond phase 
I for T> 640 K but phase IV becomes more stable at 
higher pressures. This suggests that the structure of 
phase IV needs to be revised, since theproposed Pbcn 
structure |L4| was found by Bonev et al. |4| to be thermo- 
dynamically unstable with respect to the Cmca structure 
in the whole P — T field of phase IV. 

Finally, our findings are useful in clarifying the present 
debate on the molecular character of the high P-T phases 
IV and II of carbon dioxide. As mentioned in the begin- 
ning, the C=0 bond lengths found experimentally for 
these two phases El are much larger than in phase 
I 13], which was taken as evidence that these phases 



are intermediate forms between the molecular and cova- 
lent bonded solids. Bonev et al. Q found that the bond 
lengths in either of phase II or IV structure shortened to a 
value comparable to the one in phase I; the calculated en- 
ergy difference between the long and short bond lengths 
is in excess of 3 and 6 eV per molecule respectively in 
phase II and IV. We have shown here that phase VII is a 
strictly molecular phase, with a C=0 bond length close 
to that of phase I at 1 GPa. Moreover, from the vol- 
ume and entropy jumps at the VII-IV transition, we es- 
timate an internal energy difference of about 13 meV per 
molecule between the two phases at 726 K: this strongly 
suggests that phase IV, and most likely phase II , are 
strictly molecular phases too, in agreement with a recent 
spectroscopic study Q. 



We thank B. Canny (IMPMC) and M. Mezouar (ESRF 
ID27) for their help with the experiments. 



[2 

[3. 

[4 
[5 

[6 
[T 
[8 
[9 
[10 

[11 

[12 

[13 

[14 

[is; 

[16 

[17; 

[18 



V. Iota, C. S. Yoo, and H. Cynn, Science 283, 1510 
(1999). 

V. Iota and C. S. Yoo, Phys. Rev. Lett. 86, 5922 (2001). 
C. S. Yoo, V. Iota, and H. Cynn, Phys. Rev. Lett. 86, 
444 (2001). 

S. A. Bonev, F. Gygi, T. Ogitsu, and G. Galli, Phys. Rev. 
Lett. 91, 065501 (2003). 

C. S. Yoo, H. Cynn, F. Gygi, G. Galli, V. Iota, M. Nicol, 
S. Carlson, D. Hausermann, and C. Mailhiot, Phys. Rev. 
Lett. 83, 5527 (1999). 

F. A. Gorelli, V. M. Giordano, P. R. Salvi, and R. Bini, 
Phys. Rev. Lett. 93, 205503 (2004). 

H. Olijnyk and A. P. Jephcoat, Phys. Rev. B 57, 879 
(1998). 

V. M. Giordano, F. Datchi, and A. Dewaele, J. Chem. 
Phys. 125, 054504 (2006). 

F. Simon and G. Glatzel, A. Anorg. Allgem. Chem. 178, 
309 (1929). 

B. Kuchta and R. D. Etters, Phys. Rev. B 38, 6265 
(1988). 

K. Aoki, H. Yamawaki, M. Sakashita, Y. Gotoh, and 

K. Takemura, Science 263, 356 (1994). 

F. Gygi, Comput. Mater. Sci. 10, 63 (1998). 

R. T. Downs and M. S. Somayazulu, Acta Cryst. C54, 

897 (1998). 

J. H. Park, C. S. Yoo, V. Iota, H. Cynn, M. F. Nicol, and 
T. Le Bihan, Phys. Rev. B 68, 014107 (2003). 

C. S. Yoo, H. Kohlmann, H. Cynn, M. F. Nicol, V. Iota, 
and T. LeBihan, Phys. Rev. B 65, 104103 (2002). 

We used the XDS [W.J. Kabsch, Appl. Cryst. 26, 795- 
800 (1993)] and Crystals [P.W Betteridge et al. , J. Appl. 
Cryst. 36, 1487 (2003)] softwares. 

u + and v~ come from the Fermi resonance between the 
overtone of the bending 2^2 and the symmetric stretching 
mode v\. 

Ref. reports a calculated energy difference of 
20.6 meV/molecule between the molecular Cmca and 
P42/mnm (phase II) structures. 



